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We search for the radiative penguin decays B® — * 07 and B® — > 77 in a 23.6 fb -1 data sample 
collected at the T(55') resonance with the Belle detector at the KEKB e + e~ asymmetric-energy 
collider. We observe for the first time a radiative penguin decay of the B° meson in the B° -> 07 
mode and we measure B(B° -> 07) = (57±^(stat)lii(syst)) x 1(T 6 . No significant B° — > 77 signal 
is observed and we set a 90% confidence level upper limit of B(B° -» 77) < 8.7 x 10~ 6 . 

PACS numbers: 13.25.Gv, 13.25. Hw, 14.40.Gx, 14.40.Nd 



Radiative penguin decays, which produce a photon via 
a one-loop Feynman diagram, are a good tool to search 
for physics beyond the Standard Model (SM) because 
particles not yet produced in the laboratory can make 
large contributions to such loop effects. The B® — > 07 
mode is a radiative process described within the SM 
by a b — ► sj penguin diagram (Fig. [T] left); it is the 
strange counterpart of the B — > if* (892)7 decay, whose 
observation by CLEO in 1993 0] unambiguously demon- 
strated the existence of penguin processes. In the SM, 
the Bg — > 07 branching fraction has been computed 
with -30% uncertainty to be about 40 x 10 -6 @, 0|. 
The B® — > 77 mode is usually described by a penguin 
annihiliation diagram (Fig. [T] right), and its branching 
fraction has been calculated in the SM to be in the range 
(0.5-1.0) xlO" 6 Neither B° -> 7 nor£° -> 77 

has yet been observed, and the upper limits at the 90% 
confidence level (CL) on their branching fractions are, 
respectively 120 x 10~ 6 and 53 x 10" 6 @. 

A strong theoretical constraint on the B® — > 07 
branching fraction is generally assumed due to good 
agreement between SM expectations and experimental 
results for b — ► S7 rates, such as in B + — > A*(892) + 7 
and B° -> A'*(892)S decays |, 0, M, El or inclusive 
B -> X s j decays P~D, E3). The B° -> 77 decay rate is 
constrained in a similar way (l3| . though various New 
Physics (NP) scenarios such as supersymmetry with bro- 
ken i?-parity [lij , a fourth quark generation (lq or a two 
Higgs doublet model with flavor changing neutral cur- 
rents [l6| . can increase the B® — > 77 branching fraction 
by up to an order of magnitude without violating con- 
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FIG. 1: Diagrams describing the dominant processes for the 
Bs —> 07 (left) and B® — > 77 (right) decays. 



straints on the B — > X s j branching fraction. 

In this study, we use a data sample with an integrated 
luminosity (Lint) of 23.6 fb -1 that was collected with 
the Belle detector at the KEKB asymmetric-energy e + e~ 
(3.6 on 8.2 GeV) collider [T3] operating at the T(55) 
resonance (10.87 GeV). 

The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a 4-layer silicon detector 
(SVD (|l3|)i a central drift chamber (CDC), an array of 
aerogel threshold Chercnkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to detect K\ mesons and to identify muons. The detector 
is described in detail elsewhere [19(. 

The variety of hadronic events at the T(5S) resonance 
is richer than at the T(45). B + , B° and 5° mesons are 
all produced in T(55) decay. B® mesons are produced 
mainly via T(5S) B*B* decays, with subsequent B* 
low energy photon de-excitation. The bb production cross 
section at the T(55), the fraction of B^ Bg events in 
the bb events, and the fraction of B*B* events among 
B s B s events have been measured to be, respectively, 
a T(BS) = (0.302 ± 0.015) nb [IJ, f s = (19.5±^°)% [n| 
and f B , s , = (9315)% 1- The B * S B Q S and B° S B° S decay 
fractions are small and not yet measured. 

Charged tracks are reconstructed using the SVD and 
CDC detectors and are required to originate from the 
interaction point. Kaon candidates are selected from 
charged tracks with the requirement Ck/{£k + Ar) > 
0.6, where Lk (£tt) is the likelihood for a track to be 
a kaon (pion) based on the response of the ACC and 
on measurements from the CDC and TOF. For the se- 
lected kaons, the identification efficiency is about 85% 
with about 9% of pions misidentified as kaons. 

We reconstruct mesons in the decay mode — > 
K + K~ by combining oppositely charged kaons having 
an invariant mass within ±12 MeV/c 2 (~ 2.5<r) of the 
nominal mass 
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We reject photons from ir° and r\ decays to two photons 
using a likelihood based on the energy and polar angles 
of the photons in the laboratory frame and the invariant 
mass of the photon pair. To reject merged photons from 
7T° decays and neutral hadrons such as neutrons and A£, 
we require an ECL shower shape consistent with that of 
a single photon: for each cluster, the ratio of the energy 
deposited in the central 3x3 calorimeter cells to that of 
the larger 5x5 array of cells has to be greater than 0.95. 
Candidate photons are required to have a signal timing 
consistent with originating from the same event. For the 
Bg — > 77 mode, photons are selected in the barrel part 
of the ECL (33° < 6 < 128°) and we require that the 
total energy of the event be less than 12 GeV. 

B° s meson candidates are selected using the beam- 
energy-constrained mass A4c = J (E™ m ) 2 - (p™) 2 



and the energy difference AE = E^ 1 - E^ m . In 
these definitions, E™ m is the beam energy and p™ 
and Eg^f are the momentum and the energy of the 
_B" meson, with all variables being evaluated in the 
center-of-mass (CM) frame. We select B® meson can- 
didates with Mbc > 5.3 GeV/c 2 for both modes, and 
-0.4 GeV < AE < 0.4 GeV for the B° -> 7 mode and 
-0.7 GeV < AE < 0.4 GeV for the B° 77 mode. 
No events with multiple B® candidates are observed in 
either data or Monte Carlo (MC) simulation. B* mesons 
are not fully reconstructed due to the low energy of the 
photon from the B* decay. Signal candidates coming 
from B*B* S , B* S B° S and B°B° S are well-separated in M bc , 
but they overlap in AE [9(. 

The main background in both search modes is due to 
continuum events coming from light-quark pair produc- 
tion (uu 1 dd, ss and cc). Rejection of this background 
is studied and optimized using large signal MC samples 
and a continuum MC sample having about three times 
the size of the data sample. A Fisher discriminant based 
on modified Fox- Wolfram moments (SFW [2l|] ) is used 
to separate signal from continuum background. The pro- 
cess e + e~ — > qcpj is a source of high-energy photons 
with low polar angles and can thus be a background 
for radiative B decays. Therefore, for the B® — > 07 
mode, we apply a more restrictive SFW requirement 
when the candidate photon is recontructed outside the 
barrel part of the ECL. This procedure is not used for 
the B® — > 77 mode where photons are selected only in 
the barrel. For the B® — > (fry mode, the SFW require- 
ment is chosen in order to maximize a figure of merit 
defined as N sig /^N s i g + N u dsc, where A s i g and N u d sc 
are the expected number of signal events coming from 
B*B* events and continuum events, respectively. A s i g 
and N u dsc are computed in the B® — ► 07 signal window 
(M bc > 5.4 GeV/c 2 , -0.2 GeV < AE < 0.02 GeV and 
I cos#h c i| < 0.8) and are normalized to an integrated lu- 
minosity of 23.6 fb _1 assuming B(B° -> <jrf) = 40 x 10~ 6 . 
The helicity angle 6>h c i is the angle between the B® 



and the K + in the rest frame. For signal events 
cos #hci should follow a 1 — cos 2 #hei distribution, while 
for continuum events the distribution is found to be 
flat. For the B® — ► 77 mode, we optimize the SFW re- 
quirement to minimize the 90% CL upper limit on the 
branching fraction computed by the Feldman-Cousins 
method 22j. The upper limit calculation requires two 
inputs: the number of observed events (iVo bs ) and the 
expected number of background events (iVbkg)- We as- 
sume iV b 8 = iV sig + N udsc and A bkg = N u dsc- N sig 
and Nudsc are computed in the B® — > 77 signal window 
(M bc > 5.4 GeV/c 2 and -0.3 GeV < AE < 0.05 GeV) 
assuming that B(B° S -> 77) = 1.0 x 10~ 6 . 

Inclusive bb backgrounds from T(5S) decays are stud- 
ied using MC samples having about the same size as the 
data sample. Backgrounds coming from B+ or B° de- 
cays are found to lie outside of the fit region. For B® 
decays, no event is reconstructed in the B® — > 77 mode. 
The B® — > ^77(77) decay is a potential background for 
the B® — > 07 mode and is studied using a dedicated 
MC sample. Assuming that its branching fraction is the 
same as its B° counterpart B° — > 7\*(892)°?7 [ll|, we ex- 
pect to reconstruct one B® — > (ji'ni.ll) background event. 
Considering the large B® — ► 0^(77) branching fraction 
uncertainty, this background is treated as a source of 
systematic error. 

For the B° s — > 07 (B® — > 77) mode, we perform a 
three-dimensional (two-dimensional) unbinncd extended 
maximum likelihood fit to Mb c , AE and cos(?hei (Mb c 
and AE) using the probability density functions (PDF) 
described below. 

The signal PDFs for Mb c and AE are modeled sepa- 
rately for events coming from B*B*, B*B° and B°B° 
with smoothed two-dimensional histograms built from 
signal MC events. The M hc (AE) mean for the B*B* 
signal is adjusted to the B* mass (the B*-B® mass dif- 
ference) obtained from B® — ► D~tt + events reconstructed 
in the same T(55) data sample. The Mb c and AE res- 
olutions for the B° s — * 07 (B° s — > 77) signal are cor- 
rected using a control sample of B° — > A'* (892)% events 
(e + e~ — > 77 events) recorded on the T(4S') resonance. 
Statistical uncertainties contained in these corrections 
are included in the systematic uncertainty. Continuum 
background is modeled with an ARGUS function [23[ for 
Mbc and a first-order polynomial function for AE. For 
the B° s — ► 07 mode, the signal (continuum) PDF for 
cos #hci is modeled with a 1 — cos 2 6*hei (constant) func- 
tion. The B° s — ► 0^(77) background PDF is modeled 
using MC events as the product of a two-dimensional 
PDF for Mb c and AE and a one-dimensional histogram 
for cos 6>hd ■ The likelihood is defined as 



n(E^)> 



(i) 



where i runs over all events, j runs over the possible event 
categories (signals or backgrounds), Sj is the number of 
events in each category and Pj is the corresponding PDF. 
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TABLE I: Efficiencies, signal yields, branching fractions and 
significances (Sig.) obtained from the fits described in the 
text. The first uncertainty is statistical and the second sys- 
tematic. The upper limit is calculated at the 90% CL. 



Mode 


e(%) 


5 S B 




S B . 


B' 


B (1(T 6 ) 


Sig. 


07 
77 


24.7 
17.8 


-0.7±?;g 
-4.7jt» 




is: 

-7.3 


-6 

-5 

4- 2.4 
-2.0 


t-7+18+12 
J '-15-ll 

< 8.7 


5.5 



Both fits have six free fit variables: the yields for 
the B*B*, B*B° and B° S B° S signals (S B , B ,, S B . B o and 
S B o B o), the continuum background normalization and 
PDF parameters, except the ARGUS endpoint which is 
fixed to 5.435 GeV. The branching fractions (B(B° -> 
07) and B{B° S — > 77)) are determined from the B*B* 
signal yields with the relations 



J B*B* 



B(B° S -» 77) x e 77 x iV s o x f B * B * , (2) 

2?(B S ° -» 7 ) x 6(0 -> K + K~) 
x e H x AT B o x f B , B , , (3) 



where e's are the MC signal efficiencies listed in Table U 
and N B o is the number of B® mesons evaluated as N B o = 

2 x Ltat x a b T 6 (5S) x /, = (2.8l°;|) x !0 6 - 

In the B® — > 07 mode we observe 18^5 signal events 
in the B*B* region and no significant signals in the two 
other regions. These signal yields are compatible with 
f B . B , = (93i£)% 0. We measure B{B° -> 7 ) x f a x 
f B l B , = (10.3l 3 2 j±1.3) x 10~ 6 and B(B° S 7 ) = 
(57i*giJi) x 10~ 6 with a significance of 5.5tr, where the 
first uncertainty is statistical and the second is system- 
atic. Systematic uncertainties and computation of the 
significance are detailed below. The measured branch- 
ing fraction is in agreement with SM expectations 0, 0] 
and with the measurements B(B° — > i^*(892)°7) = 
(40.1 ± 2.0) x 10" 6 and B(B+ -> K *(892)+ 7 ) = (40.3 ± 



2.6) x 10~ 6 We observe no significant B^ — > 77 sig- 
nal and, including systematic uncertainties, determine a 
90% CL upper limit of B(B° a -> 77) < 8.7 x 10~ 6 . This 
limit is about six times more restrictive than the previ- 
ous one 0, though still about one order of magnitude 
larger than SM expectations IS], |q bfland still above the 
predictions of NP models [LJ, [l5|, [l6|. The results are 
summarized in Table [J and fit projections in the signal 
windows are shown in Figs. [2] and [3] 

Systematic uncertainties are listed in Table HH The 
error on the signal reconstruction efficiency is dominated 
by uncertainty on the efficiency of the SFW requirement . 
This uncertainty is evaluated by comparing efficiencies 
in data and MC using the B Q S -> D~tt + control sam- 
ple. For the B° s — > 07 mode, we take as systematic 
uncertainty the B difference between the results of the 
nominal fit and the results of a fit where the continuum 
is parametrized with a second-order polynomial function 
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FIG. 2: Mbc, AE and cos#hei projections together with fit 
results for the B® — > <f>^ mode. The points with error bars 
represent data, the thick solid curves are the fit functions, 
the thin solid curves are the signal functions, and the dashed 
curves show the continuum contribution. On the Mbc figure, 
signals from B^B®, B^B® and B^Bl, appear from left to right. 
On the AE and cos (9i lc i figures, due to the requirement Mb c > 
5.4 GeV/c 2 only the B*B* signal contributes. The bottom 
right figure shows AE versus Mb c for selected data events. 
The dashed lines show the signal window. 




5.3 5.32 5.34 5.36 5.38 5.4 5.42 5.44 

M bc (GeV/c 2 ) 



FIG. 3: Mbc and AE projections together with fit results 
for the Bg — > 77 mode. The points with error bars represent 
data, the thick solid curves are the fit functions, the thin solid 
curves are the signal functions, and the dashed curves show 
the continuum contribution. On the Mbc figure, signals from 
B®Bs, B*B® and B*B* appear from left to right. On the 
AE figure, due to the requirement Mbc > 5.4 GeV/c 2 only 
the B*B* signal contributes. 



for AE. For the B® — > 77 mode, the limit obtained with 
the nominal continuum parametrization is found to be 
conservative. For the — > 07 mode, systematic un- 
certainties on B are evaluated by repeating the fit with 
each parameter successively varied by plus or minus one 
standard deviation around its central value. The pos- 
itive and negative uncertainty in B are obtained from 
the quadratic sum of the corresponding deviations from 
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TABLE II: Systematic uncertainties. 



Source 



b° -» ch b° s 



77 



Photon reconstruction efficiency 

Tracking efficiency 

Kaon identification efficiency 

SFW requirement efficiency 

MC statistics 

e (quadratic sum) 



2.2% 2 x 2.2% 
2 x 1% 
2 x 1.1% 

10% 10% 
0.8% 1.1% 
10.7% 11.0% 



Signal shape 


±2:296 ne ? 


;i- 


AE continuum shape 


nef 


;i- 


Bs backgrounds 


±°^% nej 


;i- 


B{<f> -> K+K~) 


1.2% 




Lint 


1.4% 




T(5S) 


5.0% 




fs 


+15 (W 
-12 /0 




N B o (quadratic sum) 


+ 16(W 
-13 /o 




JB*B% 


+ 7. 5m- 
-9.7 /o 




Total (quadratic sum) 


+21<y +21 
-20 /0 -19 


% 



the B value returned by the nominal fit. The signifi- 
cance of the branching fraction measurement is defined 
as sj2 (In £ max — In Cq ) , where £ max is the likelihood re- 
turned by the nominal fit and £q is the likelihood re- 
turned by the fit with B set to zero. Systematic uncer- 
tainties arc included by choosing the lowest significance 
value returned by the fits used to evaluate the systematic 
uncertainty. The B° s — > ^77(77) background is the only 
source of systematic uncertainty having a non-negligible 
effect on the significance. For the — > 77 mode, the 
90% CL limit, Sn m it, is computed by likelihood integra- 
tion, according to j Bumit C{B) dB = 0.9 x ft C(B) dB. 
Systematic uncertainties are included by convolving the 
likelihood function with Gaussian distributions for the 
parameters giving rise to systematic uncertainty. 

In summary, we observe for the first time a radiative 
penguin decay of the meson in the B® — > <f>j mode. 
We measure B(B° S -> <jyy) = (57±is(stat)±^(syst)) x 
10~ 6 , which is in agreement with both the SM predic- 
tions and with extrapolations from measured B + — > 
if*(892)+7 and B° -> K* (892)% decay branching frac- 
tions. No significant signal is observed in the B® — > 77 
mode and we set an upper limit at the 90% CL of 
B(B° -> 77) < 8.7 x 10~ 6 . This limit significantly im- 
proves on the previously reported one and is only an order 
of magnitude larger than the SM prediction, providing 
the possibility of observing this decay at a future Super 
B-factory 
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